Electronic Structure in Gapped Graphene with Coulomb Potential 
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In this paper, we numerically study the bound electron states induced by long range Coulomb 
impurity in gapped graphene and the quasi-bound states in supercritical region based on the lattice 
model. We present a detailed comparison between our numerical simulations and the prediction 
of the continuum model which is described by the Dirac equation in (2 + l)-dimensional Quantum 
Electrodynamics (QED). We also use the Fano's formalism to investigate the quasi-bound state 
development and design an accessible experiments to test the decay of the supercritical vacuum in 
the gapped graphene. 
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Introduction Graphene, a two-dimension (2D) 
hexagonal lattice of carbon atoms, exhibits special elec- 
tronic dispersion relation, in which electrons behave like 
massless relativistic Dirac fermions [ll,l3,i]. This special 
property leads to many unconventional phenomena, such 
as the Klein paradoxQ and the Vaselago lensing effect [H]- 
Due to its large "fine structural constant", graphene 
not only provides an exciting platform to validate some 
predictions of Quantum Electrodynamics (QED) in the 
strong field, but also provides an interesting "strong cou- 
pling" 2+1 QED model[6|. 

Recently, it was found that a substrate induced po- 
tential can break the chiral symmetries of massless 
Dirac equation and generate gaps in graphene electron 
spectrum[X]. The motion of electrons, therefore, can 
be described by massive 2D Dirac equations. It was 
predicted that a Coulomb charged impurity in gapped 
graphene behaves the same as the heavy-atom in the 
QED theory. The charged impurities have a significant 
effect on the electronic structures of the gapped graphene 
B 10, EH, El- Depending on the charge Z, the bound 
states can be induced inside the gap, and the quasi-bound 
states can also be generated when the charge Z is suffi- 
ciently larger than the supercritical number Zc- If the 
energy of bound state is above the middle of the gap, 
the wave function of bound electron state in the pure 
Coulomb potential can be well determined by using the 
continuum model. If the energy of bound state is below 
the middle of the gap, the wave function of bound state, 
however, cannot be well described because of the singu- 
lar nature of the pure Coulomb potential near r ~ 0. 
To solve this problem, a rough but simple approximation 
is adopted to remove this singularity condition by taking 
into account the finite spatial extension of nuclear charge 
with radius R based on 3+1 QED model [l^. When a 



continuum model is used to describe the behavior of elec- 
trons in graphene near the Dirac point, in fact, how to 
choose the boundary condition r ^ is still an open 
question. For example, Ref.[14] chose a "Zigzag edge" 
boundary condition to discuss the vacuum polarization 
in gapless graphene, but Ref.[lH| used an "infinite mass" 
boundary condition. This situation, however, dose not 
exist in the tight-binding model due to its lattice scale 
cutoff. 

In this paper, based on the lattice model, we use the 
tight-binding approach to numerically calculate the elec- 
tronic structure in gapped graphene. From the results of 
local density of states (LDOS), we discover that the en- 
ergy position of bound states can be determined within 
the non-supercritical region, and supercritical charge Z^ 
can also be determined using the continuum model and 
the lattice model. We also discuss the relationship be- 
tween the gap and critical charge Zc- Finally, we use 
the description of resonances to explain the quasi-bound 
state development in a supercritical regime. 

The Lattice Model and Continuum Model In 
what follows, we consider a single attractive Coulomb 
impurity placed at the center of a honeycomb lattice in 
a graphene sheet. The corresponding Hamiltonian in the 
tight-binding form is: 



(1) 
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where t = 2.7eV is hopping energy between the near- 
est neighbor atoms. The operator a^a) and b^{b) denote 
creation (annihilation) an electron on the sublattice A 
and sublattice B, respectively. The first term of Hamilto- 
nian describes the hopping between the nearest neighbor 
atoms. The second term is the on-site energy of sublat- 
tice A and sublattice B. M corresponds to the mass (or 
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FIG, 1: (Color online) Subfigures la, lb and Ic describe the 
bound states of the nearest neighbor with different gaps. Id 
shows the relationship between the energy of bound states 
and the gap when Za = 0,4, The line is the theoretical pre- 
diction for ZoL = 0,4 and j = 1/2, Subfigures le,lf and Ig 
describe the bound states of the nearest neighbor with dif- 
ferent Coulomb charge value when M = 0,05t. Subfigure Ih 
shows the relationship between the energy of bound states 
and the Coulomb charge values when M = 0,05t, The black 
line is the theoretical prediction. 
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FIG, 2: (Color online)2a, LDOS of the nearest neighbor 
atoms with different Za) when the gap M is set to be 0.05t, 
2b, LDOS at the energy scale M < E < t differs obviously 
when the bound state reaches OeV. 



states within the gap \E\ < M, as a result, satisfies the 
following condition: 



the gap) of Dirac fermions. It is weU known that the 
non-zero mass term can open a 2M wide gap within the 
band spectrum. Many conditions can lead to a band-gap 
within graphene. For instance, Ref 7] reports that the 
SiC substrate can open a gap of ~ 0.26eV in a single- 
layer graphene, which is consistent with the calculation 
of the first principle [l^. In our calculations, we select 
M from 0,05i (or O.lSeV) to O.lOt (or 0.27eV). Z reflects 
the impurity charge strength and e is the effective dielec- 
tric constant. When energy level is close to the Dirac 
point, the Hamiltonian in gapped graphene can be ap- 
proximately described by a continuum model: 

H = -ivFicrA + (J2dy - ^) + Ma3 (2) 

where ui 2,3 are the Pauli matrices. We can view a = 
as "flne structural constant" in graphene, vp — ^ta 
is the Fermi velocity of graphene and a = 1,42 A, AI is the 
mass of the Dirac fermion. Since vp is sufficiently small 
compared with the velocity of light in the QED theory, 
a ~ 1 then becomes much larger than that in the QED 
theory. The large value of "fine structural constant" of 
graphene means that the perturbation expansion, which 
works well in QED, can not directly apply to discuss the 
many body problems in graphene , 

Bound states above the middle gap In the con- 
tinuum model, the eigenfunction of a pure Coulomb po- 
tential can be described by the confluent hypergeometric 
function The regularity at r — > and r —i- 00 

requires the confluent hypergeometric functions reduce 
to polynomials at the same time. Energy of the bound 



E^, = ^^M^^ (3) 

V ^ (n+'rV 

with the root 7 = — (Za)^. Here n and j takes an 
integer value and half of an integer, respectively. That 
is, n = 0, 1, 2., and j = 1/2, 3/2.,, j is the isospin-orbital 
momentum number [iTj , Non-supercritical region is de- 
fined because 7 is always real for all angular momentum 
channels. If the control charge exceeds the value Za > j, 
the solutions of the continuum model break down since 
the root 7 becomes imaginary. Hence, all bound states 
below the middle of gap do not exist. This artifact can 
be remedied by removing the singular behavior of the 
pure Coulomb potential. When the energy of the bound 
state above the middle of gap, one can see from Eq, (3) 
that the energy of the bound state is linearly proportion 
with the gap width, and the lowest bound state reaches 
the middle of the gap E ~ OeV at a critical value that 
satisfies Za = j = 1/2, 

As shown in Fig,l, the LDOS spectrum for the bound 
state above the middle gap is calculated numerically by 
the tight-binding approach. By fixing the electromag- 
netic coupling Za = 0,4, we can plot the LDOS spec- 
trum of the bound states with different gap widthes (M) 
as shown in Fig, la to Fig, Ic, Except that the low- 
est bound states {n = and j = 1/2) in the gap can 
be clearly resolved based on the LDOS spectrum, the 
other bound states are very close to the edge of the pos- 
itive energy continuum and thus one cannot clearly dis- 
tinguish them from Fig, la to Fig,lc, The energy of the 
lowest bound state is plotted in Fig, Id as a function of 
the gap width. The dotted line is depicted based on our 
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numerical simulations in the lattice model. The theoret- 
ical results is also plotted by using the equation (3) of 
the continuum model. When the gap width is not too 
large, the energy of the lowest bound state calculated by 
the numerical simulation show a slightly deviation away 
with that predicted by the continuum model. As the gap 
M increases, the deviation becomes larger. 

If the gap width is set to be M = O.OSt, Figl.e to 
Fig.lg illustrate the LDOS spectrum of the bound states 
with different chosen electromagnetic coupling factor, 
Za. The peak of LDOS corresponding to the lowest 
bound state gets close to the middle gap as Za increases. 
Fig. Ih plots the energy of the lowest bound state as a 
function of Za. The dotted line gives the result of nu- 
merical simulations and the fitting curve is plotted using 
equation (3) of the continuum model. The numerical sim- 
ulations fit Eq. (3) well when Za is small. Fig. Ih also 
imply that the electromagnetic coupling Za needs to be 
larger than Za — 1/2 in order to make the energy of the 
lowest bound state to reach the middle of the gap. 

The effect of Coulomb charge on energy continuum is 
shown in Fig. 2. Except for the existent of bound states 
in the gap, the main feathers of the LDOS spectrum on 
gapped graphene is almost the same as that in the gaplcss 
graphene [15]. The strong re- normalization of the van 
Hove singularities are also observed, [l^ . Note that the 
LDOS spectrum near the edge of positive energy contin- 
uum appears some interesting physical signature, which 
is discussed in the following. 

Bound states below the middle gap The Dirac 
equation with the pure Coulomb potential can only pro- 
vide the solution of bound state above the middle of gap. 
When Za exceeds j = 1/2, 7 becomes imaginary in this 
continuum model, which does not have a regular solution. 
For the pure Coulomb potential, the continuum model 
points out that this is critical point when Za equals 1/2. 
In the lattice model, the bound state, however, can be 
gradually enter the energy region below the middle of gap 
as Za increases. Thus, this point seems not critical. In 
this paper, for a convenient comparison between the re- 
sults obtained by the lattice model and by the continuum 
model, we still define a critical value Z^a which repre- 
sents the energy of the lowest bound state just touches 
the middle of gap E = QeV . More interesting, our nu- 
merical simulations reflect that some physical signatures 
are associated to this point. 

Fig. 3b ~ Fig. 3d show the situation that the bound 
states enter below the middle of gap when we fix the gap 
width M = 0.05t. The bound states exist below the mid- 
dle gap as the Za increases. The relationship between 
the critical charge Z^a and the gap width M is numeri- 
cally calculated and the results are shown in Fig. 3a. The 
trigonal dots stands for the numerical results of critical 
charge Z^a. Therefore, the value Z^a depends on the 
gap width M in our lattice model. For example, the 
critical value Z^a is about 0.6 when the gap width is 
M = 0.05t. It is more interesting to note that there 
still exist some physical signatures related to this critical 
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FIG. 3: (Color online) 3a. Relationship between critical charge 
Zoct and Dirac mass. Subfigures 3b, 3c and 3d describe the 
bound states for Za = 0.55, 0.60 and 0.65 in the gap, respec- 
tively. 



value. As shown in Fig. 2b, the LDOS spectrums of the 
positive energy continuum near the gap appear signifi- 
cantly different for Za < Z^a versus Za > Z^a. When 
Za < Zga, the LDOS of positive energy continuum de- 
creases as E/t gets close to the gap region. However, for 
Za = Zga, the LDOS of positive energy continmmr near 
the gap is almost independent of E/t. When Za > Zga, 
the LDOS spectrum near the gap appears an increasing 
trend as E/t decreases. The physics behind this phe- 
nomenon needs to be explained in our future research. 

Supercritical Region and the Decay of Vacuum 
As Za continues to increase, the energy of the low- 
est bound state becomes close to the edge of band gap 
{E = —M). When Za is sufficiently larger than the criti- 
cal value Zctt, the energy of bound state will drop into the 
negative continuum spectrum. In this case, the bound 
state changes its characteristics and becomes a quasi- 
bound state. As predicted by the QED theory 13], the 
neutral vacuum changes into the charged vacuum. The 
processes with this change is called the decay of vacuum. 
We first calculate the supercritical charge Zca in gapped 
graphene by using the continuum model and the lattice 
model, respectively. To determine the diving point by 
the continuum model, one needs to remove the singular 
behavior of the pure Coulomb potential. A simple choice 
is that the potential takes the form V — < R, 

and V — — otherwise. Hence, the supercritical 
charge Zca for the lowest bound state can be determined 
through the boundary condition at r = -fifisl. [isl]: 



Ji{Zca) _ 1 
Jo{Zca) 2Zca 



(4) 



where pc = y/MlRZ^ and ly = 2^(Zca)2 - j2. This 
is a transcendental equation for Za and the numerical 
result gives Z^a = 0.78 at MR = 0.02. If the gap width 
is set to be M — O.Q5t, the corresponding R equals about 
0.6a. 

Based on our lattice model. Fig. 4a shows the munerical 
results when the lowest bound state enters into the nega- 
tive energy continuum. The results indicate that the su- 
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FIG. 4: (Color online)4a. The evolving bound states near 
E = —M when Za < Zca. We can observe that LDOS on 
the bound states becomes larger. 4b. The evolving resonances 
near E = —M when Za > Z^a. The arrow points towards 
the edge of gap E — —M. The inset subfigure 4c shows the 
relationship between ZcCi and gap width M. 

percritical charge is about Zca — 0.756 when M — 0.05t. 
Inset in fig. 4a reflects the relationship between Zca and 
the gap width M. The square dots represents the nu- 
merical results based on the lattice model. The fitting 
line implies that Zca depends linearly on the gap width 
M, which satisfies the result of the above transcendental 
equation (4), when MR « 1. 

If Za > Zca, the lowest bound state enters the neg- 
ative energy continuum and becomes a resonant state. 
The LDOS spectrum clearly shows the resonant state in 
negative energy continuum (refer to Fig. 4b). The LDOS 
spectrum also shows that the second lowest bound state 
is distinguishable. From the numerical results, we can see 
that the width of quasi-bound state becomes large as Za 
increases. This behavior can be explained by the method 
developed by U. Fano for the auto-ionization of excited 
states in Atomic Physics. We assume that new contin- 
uum wave function in the supercritical region can be ex- 
panded as \XE >= a{E)\ct)o > +J dE'bE'{E)\^PE' 
where {(po > and \tpE > represent the discrete eigenstate 
and the continuous spectrum when Za = Zca. If Za 
exceeds Zca, V = {Za — Zca)/r — 6za/r can be con- 
sidered as a perturbation potential. On the boundary of 
the supercritical region, we can simply revise the results 
in Ref.dl: 



where AE = Sza < 0o|l/?'|0o > and ry = 2tt\Ve\'^ = 
2TT{Szaf < ^E\l/r\(j)o >. The LDOS is defined as 
N{e,r) = EbI < r\xE > \^S{e - E) '14]. We obtain 
that N{e,r) relates to 100 > as |a(i?)p| < r|0o > 
Er = Eq + AE depends linearly on additional Coulomb 
charge while width of the quasi-bound rj quadratically 
increases with the additional Coulomb charge. These ob- 
servations is consistent with our numerical results quali- 
tively as shown in fig. 4b. 

Particularly important, some gedanken experiments 
were proposed to test the interesting physical process re- 
lated to the supercritical vacuum. '13] For example, the 
feature of the spectrum of the emitted positron depend- 
ing on the " critical duration" which reflects the decay of 
the supercritical vacuum. Unfortunately, there still exist 
difficulty to perform this experiment in QED. Since "the 
Compton wavelength Xg = hvp/M^^ in gapped graphcne 
is sufficiently large comparing with that in QED, one can 
use the gating technique to design experiments to mea- 
sure the decay of the supercritical vacuum in a gapped 
graphene. 

Summary In this paper, considering the contribution 
of the long-range Coulomb impurity in gapped graphene, 
we numerically calculate the LDOS spectrum using the 
lattice model. The relationship between the critical value 
Za and gap M is investigated, which corresponds to the 
lowest bound state at the middle gap or at the diving 
point, respectively. We also compare the numerical simu- 
lation results of the lattice model with that of the contin- 
uum model. The width of quasi-bound state is explained 
by the Fano's formulism. 

Noted added. While preparing this manuscript, we 
were aware of preprint [20j that has some overlapping with 
this paper. 
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